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We observe a strong correlation between magnetization relaxation and electrical resistivity in thin 
Permalloy (NigiFeig, "Py") films. Electron scattering rates in the films were affected by varying 
film thickness and deposition conditions. This shows that the magnetization relaxation mechanism 
is analogous to "bulk" relaxation, where phonon scattering in bulk is replaced by surface and defect 
scattering in thin films. Another interesting finding is the increased magnetization damping with Pt 
layers adjacent to the Py films. This is attributed to the strong spin-orbit coupling in Pt, resulting 
in spin-flip scattering of electrons that enter from the Py. 



I. INTRODUCTION 

The Gilbert form of the Landau-Lifshitz equation de- 
scribes the small angle precession of magnetization in a 
ferromagnet, 
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Here M is magnetization and 7 = g \e\ /2mc is the gyro- 
magnetic ratio, and a is the Gilbert damping coefficient 
that affects the magnitude of the viscous damping term. 
H c h is the effective magnetic field seen by the magne- 
tization, and is expressed in terms of the free energy as 
H c ff = — \7 M T. The Gilbert damping coefficient a, con- 
trols how rapidly the magnetization equilibrates in the 
absence of external stimulus. This obviously makes a a 
key parameter in the description of high speed dynam- 
ics in magnetic materials. A few areas where a plays a 
vital role are in devices that rely on faat magnetization 
reversal (e.g. in Giant Magnetoresistive-El, Magnetic Tun- 
nel Junction-a or other spintronic devices), current in- 
duced magnetization reversals, generation of microwaves 
by spin currentsa etc. 

By the 1970s it had been shown that intrinsic mag- 
netization relaxation in transition metal ferromagnets 
could be explained by electron scattering by phonons and 
magnons. The former process, mediated by thepSgku 
orbit interaction, occurs both withEn3, and withoutLltHjlJ 
the accompaniment of a spin-flip. In the former case 
a ~ 7- , where r _1 is the electron scattering rate. In the 
latter case the angular momentum relaxes as scattered 
electrons repopulate the magnetization-direction depen- 
dent Fermi volume, and a ~ r is expected. Magnon 
modes can relax through exchange interaction with a 
conduction electron, causing its spin to flip (this can be 
viewed simplistically as s-d-ejccJiange accompanied by a 
spin-flip of the s-electron)EjEjE3. The conduction elec- 
tron spin then relaxes to the lattice through the spin- 
orbit interaction. This also results in a ~ t" 1 to leading 
order. 

More recently, new effects were predicted for ultrathin 
films and for multilayers that would contribute to the 
effective a in Eq. (|j). Arias et al. showed that under 



certain conditions in ultrathin films the uniform mode 
(k = 0), excited by ferromagnetic resonance (FMR), can 
be scattered into k magnons by anisotropic sur- 
face defects, where k is the wavevectorta This two- 
magnon scattering theory was used by Azevedo et al. to 
explain experimental results of FMR linewidth and reso- 
nance field in NiFc filmsEl Bcr ger predicted that transfer 
of electron spin angular momentum between two ferro- 
magnetic layers, separated by a nonmagnetic laver, con- 
tributes to the magnetization relaxation (i.e. a)E-3. The 
experimental results of Urban et al. confirmed an increase 
in FMR-linewidth with two layers of Fe separated hv a 
nonmagnetic layer, compared to a single layer of FeQ 



We have studied FMR and electronic transport in 
NiFe-films (NigiFeig, "Py" for short), in which the elec- 
tronic scattering rates were affected over a wide range 
by: (a) changing the surface scattering contribution by 
varying the film thickness, (b) changing film deposition 
conditions, and (c) choosing different interfaces and sur- 
face treatment. We observe a strong correlation between 
the Gilbert damping coefficient and resistivity, i.e. a oc p 
in our single layer Py-samples. Our results show that the 
dominant magnetization relaxation mechanism in these 
samples involves electron scattering, and is seemingly 
insensitive to whether the scattering occurs within the 
"bulk" of the films or at the surface. This explains wiw 
a is observed to increase with decreasing film thicknesslill. 
It also implies that the effective a in magnetic devices, 
small in at least one dimension, made with transition 
metals or alloys, is expected to be considerably larger 
than that intrinsic to the bulk material, due to an in- 
creased surf ace/ volume ratio and to enhanced spin re- 
laxation at interfaces^. Our data also appear inconsis- 
tent with the two-magnon scattering theory. Further, 
we observe that with nonmagnetic (NM) Pt enclosing 
our Py- films, the magnetization relaxation is significantly 
enhanced, in addition to the electron scattering related 
mechanism above. The enhancement is attributed to 
spin relaxation of conduction electrons that leave the Py- 
layer, both at the interfaces and within the NM-layers. 
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II. EXPERIMENT 

Our NiFe films were deposited by dc-magnetron sput- 
tering in a vacuum system with a base pressure of 
2 x 10~ 8 torr. During deposition they were (all but one 
series) exposed to a uniform magnetic field of ~ 150 Oe to 
induce uniaxial in-plane anisotropy. X-ray results show 
that the Py films are (111) textured. We made two series 
of samples of Si/SiC^/Py/PR, where PR is photoresist, 
used to protect the film from oxidation. One of these 
(called o-Py) was grown in a uniform applied magnetic 
field, the other (called d-Py) without a deliberately ap- 
plied field. Both series have a clearly defined easy axis, 
the direction of which in the d-Py was presumably de- 
fined by the Earth's field. Resistivity measurements on 
these samples were made using the van der Pauw method. 
They indicate that the d-Py series has more disorder than 
the o-Py ("ordered"). We also studied the effect of ad- 
joining Py layers with 80 A thick NM-metallic layers, 
i.e. Si/Si0 2 /X/Py/X, where X is Cu, Nb, or Pt. The Py 
thickness within each sample series was varied by deposit- 
ing a terraced structure on a single wafer using a mov- 
able shadow mask. The films were then lithographically 
patterned into arrays of discs, 1 or 2 mm in diameter. 
From each deposition we thus obtained a series of samples 
of different thickness, with minimal variations in growth 
conditions. We also made one sample of Si/Si02/Py, 
1000 A thick, and ion-milled it several times, measuring 
its thickness and magnetic properties between millings. 



To obtain the Gilbert damping in our samples we mea- 
sured their in-plane magnetic susceptibility in an FMR- 
experiment with swept frequency and fixed dc magnetic 
field. The experimcntalj-setup is essentially the same as 
that of Korenivski et altB. The magnetic softness of Py 
(H c < 4 Oe) allows the experiments to be done with ap- 
plied dc-fields H < 150 Oe. Our films are thinner than 
the skin depth at the corresponding resonance frequen- 
cies, i.e. w r /(27r) < 3.5 GHz for the uniform mode of spin 
precession. The exchange stiffness in Py causes higher or- 
der spin wave modes to appear at much higher frequen- 
cies. The ac-field is considered uniform throughout the 
films and a quasistatic approximation relates the inter- 
nal and external fields. These assumptions are supported 
by the Lorentzian lineshape of our resonance peaks, and 
holds even for the thickest samples of ~ 1000 A. The 
FMR experiments were done at a small precession-cone 
angle, the ratio of the amplitudes of the ac- and dc-fields 
being ~ 10~ 4 . There was no detectable change in the 
FMR (susceptibility) when the rf-power was increased 
by 15 dB. 

The conditions above allow us to fit the susceptibil- 
ity very well with a linearized form of Eq. (|l]) . The free 
energy J- 1 includes the Zeeman energy, a demagnetiza- 
tion term, uniaxial in-plane anisotropy and a uniaxial 
out-of-plane (interface) anisotropy term. In the coordi- 
nate system shown in Fig. [j], under the assumption that 
the applied dc magnetic field H is in-plane, T can be 
expressed as, 



J 



T = -MH sin 0sin(0 + ip) + 2nM 2 cos 2 6 + K u sin 2 9 cos 2 + 2^ cos 2 
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where K u is an in-plane uniaxial anisotropy constant, 
K s = (K s i + K s2 ) /2 is a surface anisotropy constant 
representing the average anisotropy of the upper and 
lower surfaces, and d is the film thickness. The in-plane 
uniaxial anisotropy is determined by fitting the angular 
variation of the resonance frequency, u> r , in the plane of 



the samples. Typical results are shown in Fig. 0, which 
displays a fit to cj r /27r as a function of the equilibrium 
magnetization angle, 4>. 

For the free energy given in Eq. (||) the resonance fre- 
quency is given by, 
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With the applied field and the magnetization coincident 
with the easy axis of the sample, i.e. 4> = 0, ip = 9 = ir/2, 
for each series of samples we use the leading terms in the 
field dependence of u) r in Eq. (|j|) to get an initial estimate 
of 7. Similarly we use the thickness dependence of oj r to 



obtain an initial value for K s , as shown in Fig. ^| (b). 
Subsequently we obtain the Gilbert damping parameter 
a, and refined values for K s and 7 by fitting the com- 
plex susceptibility. With this approach we arrive at a 
self-consistent result where a constant accounts for 
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FIG. 1: Real and imaginary part of susceptibility as a func- 
tion of frequency of a 80Nb/477Py/80Nb sandwich struc- 
ture (numbers denote layer thickness in A) structure with 
Hdc = 60 Oe along the easy axis of the film, <f> — 0, and 
9 = ij) = 7r/2. Also shown, the coordinate system used to 
describe the free energy of the film lying in the x — y-planc 
with the easy axis in the ^-direction. 



the shift in the resonance frequency as a function of film 
thickness and where the g-value is constant for any given 
series of samples. 

Even without our careful determination of the above- 
mentioned values we could get a good estimate of a, since 
its effect on the resonance curve can not be mimicked by 
adjusting other parameters. However, as damping in- 
creases, the resonance peak gradually disappears and it 
becomes difficult to estimate a accurately, even with our 
fitting procedure. This is reflected in the size of the error- 
bars in Fig. H (c). We emphasize that in our treatment 
a is the dimensionlcss Gilbert damping coefficient and 
not the frequency linewidth. Although we observe the 
linewidth following the same trends as a as a function of 
film thickness, we believe that the fitting procedure gives 
more reliable results than measuring the linewidth. Also, 
a represents the total effective damping, including both 
intrinsic and extrinsic damping, such as surface damping. 



III. RESULTS AND DISCUSSION 

The Gilbert damping coefficients a, for the two PR- 
coated series (o-Py and d-Py) are displayed in Fig. || (a), 
as function of inverse film thickness 1/d. Values for the 
thickest films approach the bulk damping, but as the 
thickness decreases a increases dramatically. In the o- 
Py the damping doubles from the thickest film to the 
thinnest. The effect is much more pronounced in the d- 
Py, where the thinnest film has roughly six-fold the bulk 
damping value. The room temperature resistivity p(d), 
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FIG. 2: (a) The dependence of the resonance frequency on the 
in-plane magnetization angle in a 682 A thick sample from 
the o-Py series, with H^c = 90 Oe in the plane of the film. 
From a fit to the oscillatory angular dependence we obtain the 
in-plane anisotropy field, here Hk = 10 Oe. (b) Resonance 
frequency ui r versus sample thickness d for the o-Py series. 
The shift in lj 2 . scales as 1/d. This is accounted for by surface 
anisotropy, in accord with Eq. [| The line is a least squares 
fit and corresponds to K s = 0.28 erg/cm 2 . 



of these two series is shown in Fig. || (b) . The changes 
in p reflect quite accurately the corresponding changes 
in a, i.e. Ap/po ~ Aa/ao, where the subscript refers 
to thick film values. This suggests the existence of a 
simple relationship between a and p, which is addressed 
below. The reason for plotting these data as function of 
1/d is, that if one naively assumes that a can be sep- 
arated into independent bulk and surface contributions 
a 1/d-dependence is expected as the surface/volume ra- 
tio changes. This is analogous to the assumption made 
for electron scattering in the Fuchs-Sondheimer theory of 
surface scatteringc3. 

However, neither the o-Py nor the d-Py series ex- 
hibit a perfectly linear relationship between a, or p, and 
1/d, although the deviation of the o-Py samples from 
a straight line appears small. Neglecting the small de- 
viation, a Fuchs-Sondheimer-type analysis of the o-Py 
resistivity results in a bulk resistivity at room temper- 
ature of pb = 24 pQcm and mean-free-path A = 96 A 
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FIG. 3: (a) Gilbert damping coefficient a, and (b) resistivity 
p of the d-Py (disordered) and o-Py (ordered) series as func- 
tion of inverse film thickness at room temperature. The 
two series exhibit quite different thickness dependence, (c) 
displays the correlation between a and p. 



gether, as in Fig. || (c), they fall on a single curve. 
Note that the thickness of the films in the overlap re- 
gion is quite different for the two series. For instance, 
the thinnest o-Py film in Fig. || (c) is 35 A thick, and it 
corresponds roughly to a d-Py film of thickness 65 A. At 
least to leading order, a simple proportionality, describ- 
able by one constant, of experimental a to p as functions 
of d is apparent in Fig. [| (c). However, examination of 
the data in Fig. || (a) and (b) reveals that even two pa- 
rameter fits of a, or p = a + bd~ n to the data would 
not be fully satisfying. These remarks suggest immedi- 
ately that a significant contribution to viscous damping 
in very thin ferromagnetic metals is connected with the 
electron scattering giving rise to resistivity in a general 
way not depending on whether the electron scattering is 
by phonons, defects, or surface irregularities. 

Theoretical estimates by Kamberskycl are particularly 
helpful in distinguishing a low-temperature damping 
term proportional to the electron scattering time r from 
a high-temperature term proportional to r , when elec- 
tron scattering by phonons is prevalent. Electron scat- 
tering in our samples is dominated by surface and defect 
scattering. It is mainly caused by fixed electrostatic po- 
tentials associated with compositional interdiffusion and 
structural irregularities. It is therefore more appropriate 
to use the results of Heinrich et a/Jl3 to account for our 
data. They considered the effect of s-electron spin relax- 
ation on s-d exchange and magnetization relaxation: 

By definition, the equation a = A/7M relates the di- 
mensionless Gilbert damping coefficient to the Landau- 
Lifshitz parameter A. Combination of Eq. (21) in Ref. [l(] 
and the assumption that the spin-relaxation rate t~ ~ 
(V , where £ is a constant, and r _1 is the ordinary elec- 
tron scattering rate (characteristic of electrical resistivity, 
i.e. it includes both scattering events that are accompa- 
nied by a spin-flip and those that are not), yields the 
estimate, 
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(the corresponding low temperature, or residual, values 
are pb,rcs = 14 /ificm and A res = 215 A). It is clear from 
the significant departure of p, in the d-Py series, from 
linear 1/d-dependence, that the simple assumptions of 
the Fuchs-Sondhcimer model do not hold there. Addi- 
tional scattering mechanisms such as impurity and-grain 
boundary scattering should be taken into account! 2 !]. In 
any case, it is evident from the small resistivity ratios in 
both series, /C»295K/pros ranging from 1.85 for the thickest 
films to 1.1 for the thinnest films), that scattering as- 
sociated with film surfaces and defects (including grain 
boundaries) accounts for the major share of the total re- 
sistivity in all of these films. The strongest T-dependent 
contribution to p(T) at low T is proportional to T 2 , which 
is attributed to electron-electron and electron-magnon 
scatteringE 1 ]. 

When the data in Figs. (a) and (b) are plotted to- 



Here 7 is the free-electron gyromagnetic ratio, m* the 
effective mass of the s-electron, &f the Fermi wavevector. 
We have neglected a term H/ (4irM) ~ 10 -2 compared 
to unity. We now eliminate the scattering rate with the 
Drude conductivity formula a = p^ 1 = ne 2 T/m*. Fur- 
ther, we crudely estimate M ~ pL^n with p,Q = h-f/2 the 
Bohr magneton, and n the atomic density, to obtain, 



a 



3.7 x 10 5 (fkn) 



(5) 



estimated for a 3d metal or alloy. 

The solid line in Fig. || (c), obtained by a linear least 
squares fit constrained to go through the origin, has a 
slope of 3 x 10 4 (J7m). To satisfy Eq. (||) we must then 
have £ ~ 10 — 1 , i.e. r s _1 ~ lCPV -1 . For comparison 
we observe that the diffuse scattering of an electron is 
represented by a random walk process, where the spin 
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diffusion length Z s f = A (ts/t) 1 ^ 2 , A is the mean free path 
of the electron, and t s /t is the number of spin-preserving 
scattering events before the spin-is flipped. With pub- 
lished values for Py of fcf,= 55 An 2 ], and mean free paths 
Aj < 6 A and Aj = 46 Aa, we estimate the ratio {t s /t)~ 1 
to be 0.01 and 0.7 for the down and up spin bands, re- 
spectively. Our result lies well within that range. 

Alternative to the electron scattering mechanism, one 
may also consider that of two-magnon scattering^. Since 
both mechanisms involve interfacial effects, they both 
predict increased damping with decreasing film thickness. 
But, to account for the plot in Fig. || (c) one must show 
how two-magnon scattering predicts the equal damping 
in two films having unequal thicknesses yet exhibiting 
equal resistivities because of a compensating difference 
in their preparation. The surface roughness of our sam- 
ples, determined by atomic force microscopy (AFM), was 
the same 6 ± 1 A rms , both within and between the two 
series. Furthermore, the measured surface anisotropics 
of the o-Py and the d-Py samples are constant and iden- 
tical, i.e. K s = 0.28 erg/cm 2 in both cases. It therefore 
appears improbable that (a) the increase in a with de- 
creasing film thickness is caused by variations in surface 
anisotropy, and (b) the similarity with the changes in a 
and p is simply a coincidence. Also the theory predicts a 
field dependence of a, i.e. a decrease as the applied field 
is reduced since the two-magnon contribution vanishes in 
the limit of zero fielctl We observe a constant a in the 
field range 150 down to 30 Oe, at which point a starts to 
increase due to incomplete saturation of the films. 

The effect of adjoining Py-layers with different NM- 
metallic layers on the magnetization damping is shown in 
Fig. [|, with the o-Py series as a reference. For thick sam- 
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FIG. 4: Gilbert damping coefficient a, of Py adjoined with 
Cu, Nb, and Pt. As a reference we plot the o-Py series. Also 
shown is a Py-film that was ion milled from one side. 

pies the damping parameter remains constant, equal to 
the bulk value. However, as the film thickness decreases 



a increases rapidly, the rate of increase depending on the 
type of interface. In the Cu-, Nb- and o-Py series the 
damping at dp y ~ 35 A is approximately twice the bulk 
value in those samples. It is apparent that the effect on 
a of adjoining the Py films with Cu or Nb is the same as 
protecting the Py with photoresist, i.e. the fact that these 
thin metallic films meet with Py is insignificant in these 
cases. In contrast a very pronounced effect was observed 
with Pt-coated surfaces. Assuming a linear dependence 
of a on 1 jd the slope of the Pt-coated material is more 
than 4.5 times the slope for the other overlayers, and the 
a-value at dp y — 57 A is almost 4 times the bulk value. 
This effect was observed independently by Mizukami et 
alx-i. We have also found that successive thinning of a 
Py-film by ion milling led to greatly enhanced damping. 
As can be seen in Fig. | the effect of ion milling one side 
of the sample is at least equivalent to that of having two 
Py/Pt interfaces. 

The great increase in a in the ion-milled sample most 
likely arises from increased electron scattering at the sur- 
face due to surface damage caused by the ion-milling 
process. Again we observe a constant a as a function 
of applied field, arguing against two-magnon scattering 
effects similar to those observed by LeCraw et alrB in 
single crystal garnets. The Cu-coated samples had a sur- 
face roughness of (25 ± 3 A rms ). The Pt-coated samples 
were a factor of 3-4 smoother (7 ± 2 A rms ). Despite hav- 
ing much rougher interfaces than both the o-Py and the 
Pt-coated samples, the Cu-coated samples did not show 
enhanced damping over the uncoated o-Py and much less 
damping than the Pt-coated samples, as can be seen in 
Fig. || These observations confirm that the increased a 
in the Pt-coated series is not caused by surface roughness 
effects and point towards either an interface scattering ef- 
fect or to some-intrinsic property of the capping layers, or 
both. Bergero has predicted, for nonmagnetic layers in 
contact with a magnetic layer, a contribution to a from 
exchange coupling between localized magnetic spins and 
conduction electrons accompanied by spin flip scatter- 
ing both at interfaces and from interaction with phonons 
in the nonmagnetic layers through the spin-orbit inter- 
action. The latter would imply that layers with strong 
spin-orbit interaction, such as Pt, would provide a more 
effective damping than e.g. Cu, in qualitative agreement 
with, our results. Recently, a paper by Tserkovnyak et 
al.tll caught our attention. They presented a model of 
this same system based on the idea that the magnetiza- 
tion precession in the FM layer drives a spin current into 
the NM layer, where any spin imbalance is assumed to 
relax. Thus the same applies in their case, that Pt should 
provide a large enhancement due to its strong spin-orbit 
coupling and effective spin relaxation. 

We did not study the resistivity of the trilayer films 
coated with Cu, Nb and Pt, as the situation is obvi- 
ously much more complicated than in the single layer 
case. The resistivity of the metallic capping layers is in 
all our cases lower than that of Py. Separating the resis- 
tive contributions of capping layers, Py, and interfaces is 
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very difficult, making a comparison between a and p in 
this case less meaningful. Nonetheless, these results with 
adjoining nonmagnetic layers reenforce the above conclu- 
sion about the primacy of electron scattering mechanism 
in the viscous damping of very thin ferromagnetic films. 

The effectiveness of ion-milling in increasing a com- 
pared with the lack thereof in the Cu-coated samples sug- 
gests that although the Cu-coated samples are rougher 
than the very smooth PR- and Pt-coated samples, the 
roughness is insignificant as far as magnetization damp- 
ing is concerned. One may expect a much more complex 
surface after ion-milling, caused by a mixture of redeposi- 
tion of Py, surface oxidation and structural defects none 
of which are present in the as-deposited films. 

In conclusion we have confirmed that magnetization 
relaxation in ultrathin Py films is governed by processes 
that involve ordinary electron scattering. This is analo- 



gous to bulk relaxation except that the electron scatter- 
ing in ultrathin films is increasingly caused by surfaces 
and defects as the films become thinner, whereas it is 
to a greater extent caused by phonons in bulk materials. 
We have also observed increased damping in trilayers of 
Pt/Py/Pt which we attribute to strong spin-orbit cou- 
pling in the Pt-layers. From a practical viewpoint our 
results highlight the important connection between elec- 
tron scattering and magnetization damping. That should 
prove important when designing magnetic devices with a 
desired "optimal" dynamic response. 
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supported in part by National Science Foundation Grants 
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